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Dimer decimation scheme is introduced in order to study 
the kicked quantum systems exhibiting localization transi- 
tion. The tight-binding representation of the model is mapped 
to a vectorized dimer where an asymptotic dissociation of 
the dimer is shown to correspond to the vanishing of the 
transmission coefficient thru the system. The method un- 
veils an intricate nesting of extended and localized phases in 
two-dimensional parameter space. In addition to computing 
transport characteristics with extremely high precision, the 
renormalization tools also provide a new method to compute 
quasienergy spectrum. 

PACS numbers: 72.15. Rn+72.15-v 



For the Hamiltonian (|l|) that is periodic in time, the 
matrix elements of the evolution operator U, in the 
angular-momentum basis \m) with eigenvalues p = mh, 
read 

U m ,m, = exp[-2iLcos(mh)}{-i) m - m 'J„^ m ,(2K) (3) 

Here J r is a Bessel function of order r and K = K/ (2h) 
and L = L/(2h). The spectral problem U|V>w) — 
&~ \il>J), (lo being the quasienergy) involving unitary 
matrix of the above form was transformed W to a real 
lattice tight binding model (tbm) with angular momen- 
tum quantum number representing the lattice index 
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Kicked Harper model has emerged as an important 
model in the recent literature on quantum chaos 
The model exhibits both ballistic/extended as well as 
localized states in the regime where the corresponding 
classical system is chaotic. This challenges the con- 
cept of dynamical localization in non-integrable systems 
(such as kicked rotor), interpreted as the suppression of 
quantum diffusion due to classical chaos. Localization- 
delocalization transition of the model although has its 
roots in the corresponding un-kicked case is very com- 
plex with mixed spectrum and nontrivial multi-fractal 
characteristics [|]]|]. We use renormalization group (RG) 
techniques to study this richness and complexity. Our 
methodology has clear conceptual and numerical advan- 
tages over earlier methods for spotting critical states and 
locating localization transitions. 

The kicked Harper model (jjj is described by the time- 
dependent Hamiltonian 



H(t) — L cos(p) + K cos(q) ]T S(t-k). 



(1) 
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Here q,p is a canonically conjugate pair of variables on a 
cylinder. In the absence of kicking, the quantum system 
H = Lcos(p) + K cos(q), can be written as a nearest- 
neighbor (nn) lattice model. We take p — mh, and e %q as 
a translational operator for p, obtaining the well known 
Harper equation, Q 



\K(ip m+1 + ip m -i) + L cos(hm)ip m = eipj, 



(2) 



If we choose h — 2ixa, where a is a Roth number, the 
model exhibits localization-delocalization transition at 
K = L. H This transition has been analyzed by vari- 
ous RG schemes. [|| |) 
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where the coefficients £?™ are 

B" 1 = J r (K) sin[Lcos(mh) -nr/2- uj/2}. (5) 

For finite K, the tbm (^) effectively contributes only few 
terms as Bessel's function exhibit fast decay when \r\ > 
\K\. Therefore, the tbm describes a lattice model with a 
finite range of interaction denoted as b (b « K). In the 
limit of small K, L, lo, tbm reduces to the simple Harper 
equation (||) with e = fito Q . 

The renormalization tools have proven to be a great as- 
set in studying localization transitions in aperiodic one- 
dimensional (ID) nn lattice models. We use RG method- 
ology to investigate localization transitions in the kicked 
system (|l|) represented by the lattice model (|],||). We 
propose two independent RG schemes, which are gener- 
alization of dimer decimation scheme of the nn case g] 
which reduce the lattice problem to a (vector, or block) 
dimer, the discrete analog of the textbook example of 
the quantum barrier problem. The transport character- 
istics of the lattice model can be understood as due to 
quantum interference within the dimer. Furthermore, the 
decay of the coupling between the two sites of the dimer 
is shown to correspond to the vanishing of the transmis- 
sion probability: i.e, the localized phase is asymptoti- 
cally a broken dimer under the RG flow. As we discuss 
below, scaling analysis of the transmission properties pro- 
vide an extremely accurate tool to distinguish extended, 
critical and localized phases. The key idea underlying 
both renormalization schemes is the simultaneous dec- 
imation of the two central sites of the doubly infinite 
lattice — oo, ... — 2, —1, 0, 1, 2, ....oo, namely ±1, ±2 and 
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so on, after we have eliminated the site m = 0. Both can 
be applied to general tbm (Q) with arbitrary coefficients 

H ra- 
in the first scheme, which we will refer to as vector 
decimation, the decimation is done on a vectorized form 
of the tbm: the lattice model (||) is first transformed to 
a nn vector model where each site is associated with a b- 
dimensional vector $ m : $ m = ....Xi^ 1 ) where 

-^m = u m and for < r < b, 

r-l 

X r m = B™_ r U m + ^2 B'bWUm+r-k + Um-r+k] (6) 
fc=0 

Here we assumed, for simplicity, that the lattice model 
exhibits reflection symmetry B™ = B™ r . The lattice 
model (|ij can now be written as a nn vector tbm, 



$m-l + $m+l - V m $ m = 0, 



(7) 



where V m is a 6 x b matrix whose non-zero elements 
are given by, V m (k,k + 1) = 1, V m (k,k) = -B^_ x 5 k ,b, 
V m (k, b) = B™ +1 - B k n _ 1 and V m (b- 1, 6) = 2B™ - B™_ 2 . 
These elements have been normalized with B™. The 
dimer decimation on this vector tbm is described by the 
vectorization of the RG flow of the nn tbm [|| . We first 
eliminate the central site (m — 0) and then carry out 
the iterative process of decimating the two central sites. 
At the n th step where all sites with |m| < n have been 
eliminated, the tbm for m = ±n can be written as 



$„ +1 + G(n)$_„ - E(n)$„ = 
$_„_! + G(n)$„ - E(n)$_„ = 



(8) 
(9) 



where G(n) and E(n) are b x b matrices. With initial 
conditions G(l) = V" 1 and E(l) = V x - G(l), obtained 
by decimating the central site, the renormalized matrices 
G(n) and E(n) are given by the following RG flow, the 
matrix version of the dimer map 

G(n + 1) = [E(n)G- 1 (n)E(n) - G(n)]" 1 (10) 
E(n + 1) = V n+1 + [GfnjE^WGfn) - E(n)]" 1 

In the second scheme, which we will refer to as scalar 
decimation, we seek a renormalization scheme which pre- 
serves the banded form of the tbm (|]). We define the 
2b x 2b central matrix A with an initial value Aj^(0) '■= 



B{_ v j,k € {-&,. 



-1,1,2,. 



,b}, where B™ := 



-B™ — B™ m B® n+r / B® describes the renormalized tbm af- 
ter only the central site has been decimated from 
eqn. (|). Note that expressing u±\ from the central two 
eqs. and substituting for u±\ in the remaining eqs. one 
obtains tbm on a decimated lattice ±2, ±3 ... of the same 
band structure with a renormalized 2bx2b matrix Aj t k(l) 
in the center. The recursive scheme describing the n th 
step renormalization of the matrix A is 



i4 ilfc (n+l) = detL ijfc (n)/D(n), (11) 

(A' j+M (n) A' j+ ^_ x {n) A' j+1 (n)' 

L jtk (n) = A'_ 1M (n) A_ 1 ,_ 1 (n) A_ M (n) 

\ A' lM {n) Ai_i(n) Ai,i(n) 

where D{n) :— A_i ! _i(n) J 4i i i(n) — A-i : i(n)Ai_i(n) is 
a central 2x2 determinant, j + denotes the renormalized 
lattice site: j + := j + l,j — 1 if j > 0, < 0, respectively, 
and 

Aj, k {n); j,ke {-b, ...,-1,1,..., b} 
Bi +n A - j>bork>b 



A' jjk (n) := BiZ- 



k-j 



j < —b or k < —b 



Comparing vector and scalar decimation flows, the later 
although may lack the elegance formulation of the vector 
decimation, has various numerical advantages in addition 
to being very fast, as the only source of possible singu- 
larity, namely the central 2x2 determinant D(n), is easy 
to control, and the whole procedure is completely stable 
against oversizing the bandwidth b. However, two inde- 
pendent methods provides a unique advantage in con- 
firming many subtle features of the phase diagram that 
are discussed below. 

The important quantity that characterizes the trans- 
port properties is the effective coupling of the renormal- 
ized dimer. It is the ratio R of the off-diagonal to the di- 
agonal part of the renormalized lattice. In the vector and 
scalar decimation, the appropriate quantities are R v — 
tr(GG t )/tr(EEt) and R s = tr (A_A_ f )/tr (A+A++) 
where we name the offdiagonal and diagonal b x b blocks 
respectively as {A^) hk := Aj-k, and (A+)j t k = 
j,k 6 {1, ...,6} (assuming reflection symmetry A n j = 

A- n -l). 

To confirm that the parameter R is indeed related to 
the transmission coefficient of the model, we have done 
direct calculation of the transmission properties by solv- 
ing the scattering problem on a momentum lattice for 
a truncated kicked model. This is achieved by replac- 
ing kinetic energy term by Lcos(hm)9(l\I — \m\), 0(n > 
0) := l,6(n < 0) := 0, where the parameter M defines 
the size of the scattering region I m\ < M. This scattering 
model is equivalent to tbm (^JE)) where L is replaced by 
L6(M — |m|). Outside the scattering region, m > \M\, 
the wave-function is a superposition of properly nor- 



malized plane waves if) 



±(0 



I sin(Kj) 



-1/2 



exp(±jK;m), 



with cos(k;) — (lu + 2-kI)/{2K), I integer. The re- 
flection and transmission matrices R and T are deter- 
mined by matching the ansatze for the asymptotic so- 
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lutions, u m 

and u m — J2i> Tu'ipm^ \ for m > —M, on tbm (|]) for 
\m\ < M. Importantly, decimation scheme makes the so- 
lution of the scattering problem for large M very efficient, 
as n = M — b iterates of the RG map (11) are performed 
first in order to maximally reduce the size of the scat- 
tering region, as A(n) does not depend on truncation for 
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n < M — b. Our numerical calculations confirmed that 
results obtained using the total transmission probability 
P(M = n + b) := \T W \ 2 are fully consistent with the 
ones based on the ratios R v (ri) and R s (n). 

It turns out that the scaling exponent (3 — 
linijj-joo (3{ri), (3(n) = log S(n)/ log n (where S = P, 
R S} or R v ) provides a very effective means to describe 
transport properties. The extended states are described 
by (typically monotonic) convergence of (3{n) — > 0. In 
case of exponential localization, (3(n) — ► — oo. The de- 
cay 5*(n) <~ exp(— 2n/£) can be used to calculate the 
localization length £. In contrast, the critical states are 
characterized by negative /3 exhibiting non convergent, 
oscillatory behavior. Therefore, asymptotically, both the 
critical and the localized kicked model describes a broken 
dimer. In fig.l, we display RG flow at the Fibonacci iter- 
ates F(f), the successive denominators of the continued 
fraction approximation of a = — l)/2, which is kept 
fixed throughout the paper. However, the method can 
be implemented for arbitrary a, By iterating RG equa- 
tions upto F33 = 5,702,887, the localization thresholds 
are determined almost to machine precision. Most of our 
studies describing the variation in the transport prop- 
erties are carried out for ui = , which appears to be 
the eigenvalue for all parameters. state. However, as we 
discuss later, the RG methodology can also be used to 
compute the quasienergy spectrum and analyze its trans- 
port characteristics. 



We confirm the critical nature of the state along the 
symmetry axis K = L where the exponent (3 is always 
found to be negative and bounded. For any value of K, 
RG flow displays the oscillatory pattern shown in fig.l 
with only variation being the amplitude of the oscillation. 
Fig. 2 shows variation in (3 with K where smooth varia- 
tions are intermitantly intrupted by rapid variations. As 
we discuss below, this may be related to the existence of 
critical manifolds off the symmetry line. 
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FIG. 1. The RG flow at Fibonacci iterates, F(f). The 
monotonic curve with dots describes extended states for 
K = 4 and L = 1. The oscillatory pattern characterizing 
the critical states correspond to unkicked model (diamond) 
, kicked model with K = L = 4 (crosses). In addition, 
critical state off the symmetry line is shown for K = 4.4, 
L = 2.516711 and its dual (short and long dashed lines with 
stars) . 



FIG. 2. Variations in /3 along the line K — L in kicked 
Harper. Three RG iterates corresponding F23, F24 and F25 
(darker line) are shown. 

One of the most remarkable and intriguing result of the 
RG analysis is the cascade of localization-delocalization 
transitions off the symmetry line where the extended and 
the localized regimes display intricate pattern in two- 
dimensional parameter space as shown in fig. 3. The ex- 
tended phase of the unkicked model (K > L regime) is 
landscaped by various patches of localized regimes and 
similar behavior is seen in the localized phase consisting 
of patches of extended regime which is the dual image 
of the extended reentries. A detailed numerical analy- 
sis suggests that the critical manifolds in {K,L) space, 
corresponding to singular continuous states, can only cor- 
respond to the localization-delocalization boundary, and 
hence we do not find any fat critical phase. Further- 
more, comparing figs. 2 and 3, the rapid variations in (3 
on K = L may have its origin in collisions of the critical 
manifold with itself on the symmetry line. 

Scanning one-parameter space, say by varying param- 
eter K, one sees a cascade of reentrant transitions where 
ballistic (localized) transport reappears and survives in a 
finite window in parameter after becoming localized (ex- 
tended) as shown in fig. 4. We believe that these "transi- 
tions are related to the cascades of transitions predicted 
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by semiclassical methods ||. The reentrancies exist for 
a band of quasiperienergy states and hence will result 
in modulations in transport properties for a wave packet 
consisting of a superposition of many quasienergy states. 




1 2 3 4 5 6 

K 

FIG. 3. Two-dimensional phase diagram for uj — state 
where the shaded parts describe extended (localized) regimes 
for L > K ( K > L). Almost all intricate details have a 
dual image and this provides an additional confirmation of 
the numerics. 
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FIG. 4. Series of localization-delocalization transitions in- 
dicated by p changing from — oo to for a fixed L = 4.2 using 
F25. All transitions seen in this figure can be mapped with 
those shown in figure 3. 

Treating w as a running parameter , RG flow together 
with the duality transformation (K, L) — > (L, K) (which 
maps extended states to localized states and vice versa 
as confirmed by detailed numerics) can also facilitate an 
accurate method to determine the quasienergy spectrum 
of the kicked Harper and related models. ( see fig. 5) In 
RG scheme, it is in general difficult to distinguish lo- 
calized and forbidden values of oj as it requires knowl- 
edge of df3/duj. However, by exploiting duality, local- 
ized spectrum can be easily separated from the spectral 



gaps as the later corresponds to (3(n) -00 as well as 
/3 dual (n) -► -oo. 
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FIG. 5. Quasi-energy spectrum for fixed K = 4.2. Lighter 
dots represent extended states. Darker crosses are the lo- 
calized states which were obtained as the extended states of 
the dual model. The localization-delocalization transitions for 
various quasienergy states are clearly seen. For ui = state, 
the transitions can be mapped to those of fig. 3 and the dual 
of the fig. 4. 

In summary, RG tools uncover many complex features 
of the kicked model which would have been impossible 
with previously used methods. The dimer decimation 
approach can be applied to a variety of problems which 
include unitary models without reflection symmetry and 
can also be implemented to compute spectrally global 
quantities. We hope that these tools will provide a new 
direction in resolving various important issues underly- 
ing the frontiers of localization phenomenon in complex 
systems. 
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